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Abstract 

In order to describe excitable reaction-diffusion systems, we derive 
a two-dimensional model with a Hopf and a semilocal saddle-node ho- 
moclinic bifurcation. This model gives the theoretical framework for the 
analysis of the saddle-node homoclinic bifurcation as observed in chemical 
experiments, and for the concepts of excitability and excitability thresh- 
old. We show that if diffusion drives an extended system across the ex- 
citability threshold then, depending on the initial conditions, wave trains, 
propagating solitary pulses and propagating pulse packets can exist in 
the same extended system. The extended model shows chemical turbu- 
lence for equal diffusion coefficients and presents all the known types of 
topologically distinct activity waves observed in chemical experiments. In 
particular, the approach presented here enables to design experiments in 
order to decide between excitable systems with sharp and finite width 
thresholds. 



1 Introduction 

A system is called extended if it consists of many individual similar subsys- 
tems distributed in space, and all the subsystems have the same dynamics. 
Examples of extended systems are obtained by assembling similar subsystems 
through diffusive coupling. If all the local subsystems are in the same state, 
that is, synchronized, the extended system behaves globally as the local sub- 
system. However, if one of the subsystems is desynchronized, we can have the 
emergence of collective behavior, and the dynamic behavior of the extended sys- 
tem is dependent on both the initial spatial distribution of the desynchronized 
subsystems and the parameters of the local system. 

Here we are interested in the emergence of collective behavior of an extended 
system when the initial conditions of a local subsystem are changed. When the 
local dynamics has a stable steady state in the phase space and is close to a ho- 
moclinic orbit, we show that, in the extended system, wave trains, propagating 
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solitary pulses and propagating pulse packets can develop. These types of wave 
patterns depend on the distribution of the initial states of the extended system. 

An extended system is excitable if after being slightly locally perturbed 
from a stable steady state it returns to the quiescent state, but shows traveling 
activity waves if perturbations are large enough. Systems responding in this way 
to perturbations are said to have a threshold. Systems with excitable behavior 
appear in extended chemical reactions pQ, in the propagation of action potentials 
along nerve axons (2113) an d in population dynamics models 0J|S]. 

In the context of extended chemical reactions, there are two main types 
of models showing excitability. One of these models, the nullcline model, has 
been proposed by Tyson and Fife P and Tyson and Keener [7] to explain 
excitability effects in chemical experiments. In this model, excitability is related 
to the particular behavior of the nullclincs in the two-dimensional phase space 
of a differential equation: one nullcline is N-shaped, the other is monotonic and 
they intersect at only one point. Excitability develops as a continuous process 
when initial conditions are changed across the N-shaped nullcline, leading to an 
abrupt continuous change of the orbit length in phase space. In this case, the 
system does not have a precise threshold, but a finite narrow threshold region. 
These models of excitable systems have only one stable fixed point and are 
usually associated with two-variable models, with a slow and a fast dynamics 
8 . In the literature of excitable systems they are called Type I excitability 
models. 

Another model for excitable behavior has been proposed by Noszticzius, 
Wittmann and Stirling [Hj, McCormic, Noszticzius and Swinney ^H] in the con- 
text of the observed behavior in Bclousov-Zhabotinskii type chemical systems, 
and by Ermentrout and Rinzel in the context of propagation of action po- 
tentials along nerve axons. This model gives a qualitative explanation for the 
appearance of an excitable steady state in chemical experiments and in the 
propagation of action potentials along nerve axons through a bifurcation in the 
local kinetics, from a limit cycle to a saddle-node pair of fixed points. Some 
authors called it saddle- node infinite period or SNIPER model JU|, or Type 
II excitability model. This qualitative model has three fixed points in phase 
space, one saddle point, one stable node and one unstable focus. In this case, 
the phase portrait has a sharp threshold of excitability associated with the sta- 
ble manifold of the saddle point. The threshold of excitability has two branches, 
one connecting a saddle point with the unstable focus, and the other starting 
at the saddle point and going to infinity. Several reaction-diffusion models with 
these features appeared in the literature |12U13) . These excitability models can 
produce propagating pulses in extended media as observed in experiments |T3] . 

The SNIPER or saddle-node homoclinic (SNH) bifurcation has been pre- 
dicted by Leontovich ^H] and Andronov et al. \n5\, as a topological possibil- 
ity for the break of hyperbolicity of limit cycles in two-dimensional continuous 
dynamical systems, and is one of the simplest codimension-1 semilocal bifur- 
cations As the SNH bifurcation is nonlocal, local analysis describes the 
local saddle-node bifurcation, missing the appearance of the limit cycle |18j . 
Therefore, it is difficult to build a sufficiently general explicit analytical model 
for this global bifurcation. 

For several chemical systems in stirred tank reactor experiments, Maselko 
and Epstein |19j have constructed qualitative experimental phase portraits com- 
patible with the experimental behavior and a saddle loop bifurcation scenario. 
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In this approach, the qualitative phase portraits have one stable and two unsta- 
ble fixed points, but are obtained through a saddle loop codimension-2 (Bogdanov- 
Takens) bifurcation of a fixed point [20]. To have bounded motion in phase 
space, we must have an additional stable node external to the saddle loop. 
However, the phase portraits of these systems are topologically equivalent to 
those obtained through the saddle-node homoclinic (SNH) bifurcation. Hence, 
an excitability model with threshold can be obtained by two different bifur- 
cation scenarios, but at the excitable regime, they have the same phase space 
structure. In the following, we adopt the simplest bifurcation scenario through 
the codimension-1 SNH bifurcation. 

Our aim here is to derive a sufficiently general excitability model compatible 
with the above mentioned experimental results. We begin by constructing an 
explicit model that exhibits a Hopf and a saddle-node homoclinic bifurcation, 
both driven by independent parameters. Then, we add diffusion to the two 
phase space coordinates, obtaining a general model of a reaction-diffusion ex- 
citable media. With a bench-marked numerical method, we derive properties 
of the solutions of the one- and the two-dimensional reaction-diffusion systems. 
The Hopf bifurcation will be essentially associated with the existence of stable 
oscillations of the local kinetics, as observed in continuous stirred tank reactor 
experiments |2T]. The saddle-node homoclinic bifurcation will be associated 
with the break up of the periodic behavior of the local kinetics and the appear- 
ance of a new steady state with the excitability property 

We show that, at and beyond the SNH bifurcation, an extended chemical 
system coupled by diffusion develops traveling wave trains, propagating pulses, 
pulse packets and fronts, and aperiodic wave trains or turbulent states. Periodic 
wave trains, pulses and pulse packets are tuned by different initial conditions. 
Turbulent states and fronts are obtained by changing the parameters of the 
local model. All these situations occur for equal diffusion coefficients in the two 
phase space variables. 

In the next section, we derive the main properties of locally excitable two- 
dimensional models associated with the SNH bifurcation, we derive its geomet- 
ric features in phase space, and we characterize its bifurcation set. In section 
3, the diffusion coupling is introduced, and the properties of the solutions of 
the associated one- and two-dimensional reaction-diffusion system are analyzed 
by numerical methods. As we are dealing with a parabolic partial differential 
equation, we emphasize our analysis on the pattern formation features associ- 
ated with the initial conditions of the extended system. Finally, in section 4 we 
discuss the main conclusions of the paper. 



2 The saddle-node homoclinic (SNH) bifurca- 
tion 



In order to construct an analytical model exhibiting a Hopf and a SNH bifur- 
cation, the starting point will be the versal unfolding of the Hopf bifurcation in 
polar coordinates, |2U|: 



dr 
It 

de 
It 



[3 + br 2 



r(y + ar 2 ) 



(1) 
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where a, b , v and (3 are real parameters. In the literature of chemical kinetics, 
equation is the Stuart-Landau equation with real parameters |22j . 

If, a < Q, j3 ^ 0, and v > 0, equation (|TJ has one limit cycle (in Cartesian 
coordinates) in phase space and an unstable focus at the origin. The radius 
of the limit cycle is ro = y/—v/a. For, a < 0, (3 ^ 0, and v = 0, the system 
(2.1) has a supercritical codimension-1 Hopf bifurcation. For, a < 0, and ^ < 0, 
the origin in phase space is a stable focus. If (3 = and z/ = 0, we are in the 
conditions of a (Bogdanov-Takens) codimension-2 bifurcation [201 ■ In this case, 
Q is no longer a versal unfolding for this bifurcation. 

One of the features of the normal form JU is that the phase equation is 
phase independent. From the analysis of the vector field Q, for a < 0, and 
v > 0, in Cartesian coordinates, the rotation frequency of orbits in phase space 
near the origin approaches (3. On the limit cycle, the rotation frequency is 
uj — (3 + 6tq — (3 — bvja. If, a < 0, (3 > 0, and b > (3a/v, rotations near 
the origin and in the limit cycle are counter-clockwise. If, a < 0, (3 > 0, but 
b < f3a/v, the rotation near the limit cycle becomes clockwise. For, a < 0, 
(3 > 0, and b = (3a/ v, the rotation frequency on the limit cycle becomes zero, all 
the points on the limit cycle are fixed points, and the phase space trajectories 
on the limit cycle do not rotate. Under this condition (b = (3a /v), we have a 
line of fixed points. Due to this property of the normal form l[T|l. if we make the 
parameter b phase dependent, we can incorporate a saddle- node bifurcation into 
the limit cycle, obtaining a non-local, SNH, bifurcation for the full dynamics. 
So, we change the system of equations to the new system, 

riv + ar 2 ) 

(3+(b + af(9))r 2 (2) 

where a is a new parameter, and f(0) is a non-negative function that intro- 
duces a phase dependence on the phase equation. To maintain the continuity 
of trajectories in the two-dimensional phase space, we consider that f(9) is con- 
tinuous and non-constant, and /(0) = /(27r). To simplify the stability analysis, 
we further consider that f(9) has continuous derivatives. 

To derive the qualitative structure in phase space of the orbits of the differ- 
ential equation (0), we consider the equations for the fixed points, 

r(v + ar 2 ) = 
p+(b + af(6))r 2 = 0. (3) 

If, a < 0, v > 0, and (3 + (b + af(6))r 2 > or /3 + (b + af(9))r 2 < 0, for all 
9 £ [0, 27r] and r 6 [0, ro], then system © in Cartesian coordinates has one 
stable limit cycle and the origin is an unstable focus. As in the normal form 
case, the supercritical Hopf bifurcation occurs for v = 0, provided (3^0. 

If, a < 0, v > 0, and (3 + {b + af(9))r 2 = 0, for some X , 2 , . . . G [0, 2tt], 
with f'(0i) 0, then the system (2.2) has an even number of fixed points with 
coordinates (ro,#i), (^0,^2), ••■ ■ The Jacobian matrix of equation (0) at (ro,9i) 
is, 

J =(-2i/ 02r Q (b + af(9 l )) -<ff'(9 t )) (4) 
The eigenvalues of the Jacobian matrix J3J are X r — —2v and Xg — —oa>f'{6i)/a, 



dr 

dt 
(19 
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corresponding to the eigenvectors, 



e Xr = (avf (Ofi/a- 2v)e r +2r (b + af(p i ))e e 

e\ e = e e . (5) 

For, a < 0, and v > 0, the angular direction e\ e , taken at these fixed points, 
is tangent to the circle of radius ro, which, by is an invariant set in phase 
space. Therefore, the circle r — ro is the union of the stable and unstable 
manifolds connecting the fixed points with coordinates (ro,#i). According to 
the sign of the derivatives f'(9i), the fixed points (ro, 0i) are alternately saddle 
points and stable nodes, along the circle r = r . If f'(0i) > 0, then (ro, 0i) is a 
saddle. By 10, if {avf'{6i)/a — 2u) ^ the stable manifolds associated with the 
saddle points are tangent to the directions e\ T , being transversal to the circle 
r = ro, and make heteroclinic connections with the origin of coordinates. 

Solving the second equation in JSJl for f(6), we obtain, f(6) = — ((3 + 
b r o)/ ( ar o)i an( i there exists always a choice of a such that a new fixed point ap- 
pears at r = ro. As f(6) is continuous and non-constant in the interval (0, 27r), 
the new fixed point appears for, 

aS r o 2 max,/(0) l ° J 

or, 



P + br 
r 2 min e f{0) 



«s = - 2 1' " (7) 



where the maximum and minimum are taken around the local extremes of /. 
For example, if f(9) has only one local maximum in [0, 2tt], the new fixed point 
appears for a = as given by equation ©. If f(0) has only one minimum, the 
new fixed point appears for a = as given by equation i ff} . 

The appearance of these fixed points on the limit cycle, when we vary the 
parameter a across the value as, is the SNH bifurcation. At the SNH bifur- 
cation, we have locally in phase space one or several saddle-node bifurcations 
depending of the number of local maxima and minima of /(#). 

The SNH bifurcation is tuned by the parameter a. For example, choosing 
the function f(6) = sin 2 (#/2), and solving the equations for the fixed points, 
|J3J|, we obtain, 

asin 2 (#/2) = -4 -b=^-b 
and the SNH bifurcation occurs for, 

as - & — b (8) 
v 

where one non- hyperbolic fixed point appears and has polar coordinates (ro, 7r). 
For, a < as-, there are two new fixed points and their stability properties are 
calculated from the eigenvalues of the Jacobian matrix 

Changing coordinates from polar to Cartesian, X = rcos(0) and Y = 
rsin(0), and with f(8) = sin 2 (0/2), system © can be written as, 
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Figure 1: Phase space structure of the differential equation © at the SNH 
bifurcation (a), and after the bifurcation (b). The parameter values are: a = — 1, 
b = 1.0, v = 1, and in a) a = —2.0 and j3 = 1; b) a = —2.1 and /3 = 1. 
According to (JHJ), the SNH bifurcation occurs for as — —2.0, corresponding to 
case a). Parameter values for cases c) and d) are: c) a = —1.0 and (3 = 0; 
d) a = —1.05 and /3 = 0. The fixed point x± is a stable node and X2 is a 
saddle point. The threshold of excitability coincides with the stable manifold 
of the saddle point, which makes a hctcroclinic connection with the origin of 
coordinates. 



G 



(X 2 + Y 2 ) (aX-(b+ ~(1 - ; X ))y) 
(3X + vY 

(X 2 + Y 2 )(aY + (b + ^l- 7 =jL=))Xy (9) 

At the SNH bifurcation, system © has one fixed point away from the origin of 
coordinates, (X = —ro,Y = 0). This fixed point is non- hyperbolic and the circle 
of radius ro is a one-dimensional center manifold. The stable manifold of this 
fixed point makes a heteroclinic connection with the unstable fixed point at the 
origin of coordinates. The fixed point (X = —ro, Y — 0) is semistable, Fig.^i). 
For, a < 0, (3 ^ 0, v > 0, and a < as, system © has one unstable focus at 
the origin, and a saddle- node pair of fixed points on the circle r — ro, Fig.^ 3 ). 
The stable manifold of the saddle point makes a heteroclinic connection with 
the unstable focus at the origin of coordinates. For, a < 0, v > 0, and a > as, 
system © has a stable limit cycle in phase space. 

At the SNH bifurcation and beyond it, the shape of the heteroclinic connec- 
tion between the saddle point and the origin of coordinates can be controlled, 
independently of the radius of the invariant circle in phase space. By J5J and 
(j2Jl, the tangent space to the stable manifold of the semistable fixed point at 
a = as is generated by the vector, 

2(3 

e\ — —2ve r eg (10) 

ro 

and the angle between the tangent space to the stable manifold and the tangent 
space to the center manifold is, 

a sc = arctan . (11) 

P 

For = 0, a sc = 7r/2, and the negative part of the X-axis is invariant for the 
flow ijjj}, Fig. [TJ;), and the heteroclinic connection coincides with the interval 
[— ?"o,0] of the X-axis. In Fig. QJi), the topological structure of phase space for 
a < as and (3 — is shown. Increasing (3, the heteroclinic connection winds 
around the origin of coordinates in phase space. 

In Fig. |3 we show the bifurcation diagram for equation I©, for the codi- 
mension-1 Hopf and SNH bifurcations, with (3 > 0. 

An important property of the system of equations © is that, away from 
the Hopf {v = 0), SNH (a = a s ) and Bogdanov-Takens [v = Q, /3 = 0) sin- 
gularities, all the fixed points are hyperbolic and there are no saddle-to-saddle 
connections. For a > as, the limit cycle is also a hyperbolic closed orbit. These 
characteristics of the flow imply that, away from bifurcations, the generic mod- 
els 10) and are structurally stable in the C 1 norm. From the applied point 
of view, this is equivalent to say that this model is robust. If (3 ^ 0, the SNH 
bifurcation is one of the three possible semilocal codimension-1 bifurcations of 
two-dimensional vector fields, and is unfolded by the parameterization of the 
vector field J!|J) in the sense that when we vary the parameter a across as, on 
both sides of the bifurcation, we have structurally stable (Morse-Smale) dy- 
namical systems |TJ|, Fig. If (3 = and v = 0, we are in a scenario that 
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Figure 2: Bifurcation diagram of the Hopf and SNH bifurcations for equation 
10, for parameter values a = — 1, b = 1 and = 1. The SNH bifurcation curve 
is given by ©■ 



corresponds to a Bogdanov-Takens codimension-2 singularity of the origin, but 
equation @ do not unfold it |2Uj . 

The differential equation @ is an analytical model that makes precise the 
concepts of excitability and threshold of excitability as defined by Winfree |23|: 
"A reaction is excitable if it has a unique steady state that the system will 
approach from all the initial conditions, but there exists a locus of initial con- 
ditions near which either of two quite different paths may be taken toward the 
unique steady state. If one of these paths is a lot longer than the other, then 
the system is excitable." The threshold of excitability is defined by the stable 
manifold of the saddle point. 

3 Traveling waves, solitary pulses and chemical 
turbulence 

It is well known that when the local kinetics of a reaction-diffusion system has 
a limit cycle in phase space, the extended system can support traveling wave 
trains |2Ij. I n the special case of local kinetics of the form (2.1) (X-oj systems), 
Kopell and Howard |24| have shown that, for small wave numbers and small 
amplitude wave trains, the spatial shape of wave trains follows the temporal 
behavior of the local coordinates along the limit cycle in phase space. This 
phenomena is usually associated with the self-oscillatory characteristics of the 
local kinetics, and therefore to target patterns and pacemakers 4 . 

The simplest way to obtain numerically stable traveling wave solutions in 
reaction-diffusion systems with self-oscillatory dynamics, is to prepare the ex- 
tended system with initial conditions at the unstable steady state of the local 
system and perturb it |25| . In this traveling wave spreads from the per- 

turbation point. However, if the initial state of the extended system has values 
on the limit cycle, any sufficiently small perturbation of the extended system 
dies out and, after some transient time, the system oscillates globally. 

We now investigate the persistence of wave type phenomena when the SNH 
bifurcation is parametrically crossed. 
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Introducing diffusive coupling into model equation @ , we obtain the system 
of nonlinear parabolic partial differential equations, 

D X V 2 X + vX- (3Y 

(X 2 + Y 2 ) (aX-(b+ ~(1 - ; X ))y] 
D Y V 2 Y + PX + vY 

(X 2 + Y 2 )(aY+(b+^-(l- - X ))x] 

(12) 

where V 2 = (j^j, or V 2 = (-§^i + g^r) is the Laplacian operator, and Dx 
and Dy are diffusion coefficients. In the following we always consider the case 
D X =D Y > 0. 

To analyze the development of wave type phenomena when we cross the 
SNH bifurcation, we integrate numerically equation l)12[l. using the bench- 
marked explicit integration method developed in Dilao and Sainhas ;2J>1 1 where 
optimal convergence to the solution of the continuous system is achieved if 
(Aimax{£>x, D Y }/{Ax) 2 ) = 1/6. We consider that the reaction-diffusion 
equation (|12H is defined over a finite one-dimensional or two-dimensional circular 
domain with zero flux boundary conditions. Both domains have the dimensional 

space length, 

L = Ny/6Atma,x{D x ,D Y } = NAx (13) 

where N is the number of lattice sites along each spatial direction. To fix 
parameters, we take equal diffusion coefficients, Dx = Dy = 2 x 10~ 5 , a = — 1, 
v = 1, and 6=1. For these parameter values, by JHJl, the SNH bifurcation 
occurs for as = —j3 — 1. 

To validate the numerical method, we have made several simulations for de- 
creasing values of the integration time step At, comparing concentration profiles 
along an one-dimensional domain. For integration times At < 0.02 and during a 
finite time interval, the profile of the computed solutions coincide over the one- 
dimensional domain, within a maximum global error in the concentrations of the 
order of 10 -2 . Decreasing At, the global numerical error also decreases, insuring 
convergence. In all the simulations, we have used the time step At = 0.01, and 
Ax = 1.095 x 10~ 3 , calculated from ifT^, 

Concerning initial conditions, we choose uniform concentrations of X and Y 
along the one-dimensional and two-dimensional extended region with a pertur- 
bation at the center of the domain. Then, at and beyond the SNH bifurcation 
(a < ag), we take two particular cases. In the first case, we consider that 
the initial state of the two variables X and Y take the values of the unstable 
focus of the local kinetics along the extended system, and we introduce a local 
perturbation at the center of the extended domain. The values of the X and 
Y perturbed values are the coordinates of the stable (a < as) or semistablc 
(a — as) steady state of the local kinetics. 

In the second case, the initial state of the extended system take the values of 
the stable or semistable steady state of the local kinetics, and the local pertur- 
bation in X and Y take the values of the unstable focus of the local kinetics. For 



dX 



dY 
~dt 
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zero flux boundary conditions, if we choose a steady state of the local kinetics 
as an initial state of the extended system, this steady state is also a solution 
of the reaction-diffusion equations, and therefore, for a finite time interval, the 
presence of boundaries have no effect on the time evolution of the concentra- 
tions near the region of the local perturbation. Under these conditions, all the 
subsequent analysis will be valid for finite and infinite domains. 

The numerical integration of equation (|12|l shows that if the local system in 
l|12l) is self-oscillatory (a > as), and if initial conditions correspond to uniform 
concentrations taken at the unstable steady state, with a perturbation away 
from it, the system develops a traveling wave train propagating along space. 
However, due to zero flux boundary condition, after a long integration time, the 
whole extended system develops bulk oscillations. A similar situation has been 
reported by Auchmuty and Nicolis for the Brusselator model. 

For parameter values corresponding to the SNH bifurcation of the local dy- 
namics (a — as — —(3— 1), we have a semistable steady state with a well defined 
threshold, Fig.^,) (and^J;)). Following the same strategy as in self-oscillatory 
systems, we prepare the initial state of the extended system at the unstable fixed 
point at the origin of phase space, and perturb it locally to the semistable state. 
In this case, the solutions of the reaction-diffusion system (|12|l develop fronts, 
periodic traveling waves and aperiodic signals, depending on the parameter (3, 
Fig- EJi)-c). For (3 = 0, Fig. EJl) , the threshold of excitability coincides with the 
negative part of X axis, and the system simply develops a propagating front. 
By increasing j3, the front becomes a moving source initiating traveling 

waves. These traveling waves propagate in the opposite direction of the moving 
front, Fig. Et>)- For even larger values of f3, the waves generated by the moving 
front become aperiodic, Fig. Ot) . 

Changing the initial conditions for an initial state at the semistable fixed 
point and a perturbation of large width (160 lattice sites) with values at the 
unstable focus at the origin of phase space, the reaction-diffusion system evolves 
in time in a completely different way, Fig. G2l)-f). For (3 = 0, Fig. \3fi), the 
reaction-diffusion system evolves to the semistable state. For larger values of 
(3, the reaction-diffusion system develops a propagating pulse, as observed in 
experimental systems, Fig. Ot), and for even larger values (3 an aperiodic state 
develops, Fig. |3f)- 

The comparison between cases a)-c) and d)-f) of Fig. leads to the con- 
clusion that target patterns, as observed in two-dimensional reaction-diffusion 
systems, are not necessarily related with the oscillatory characteristics of the 
local systems. 

To analyze the properties of periodic wave trains at the SNH bifurcation, 
we have calculated numerically the wavelengths A between consecutive maxima 
along the extended media, Fig. ^ For increasing values of /3, the aperiodic 
regime is reached. The wavelength as a function of (3 behaves as A ~ a + 
d//3. At the SNH bifurcation, the stability of wave trains seems to follow the 
same criteria as in self-oscillatory systems, where, for increasing wave number 
k = 1/A, the stability of periodic wave trains is lost 0J[^], being followed by 
an aperiodic regime. At the aperiodic regime, the distribution of wavelengths 
becomes continuous. In the literature of chemical kinetics, this phenomena is 
called chemical phase turbulence 

To analyze the persistence of wave type solutions when the SNH bifurcation 
is crossed, we have fixed (3 — 0.5. In this case, the SNH bifurcation occurs 
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Figure 3: Solutions X(x,t) and Y(x,t) of the reaction-diffusion equation l|12JI 
at the SNH bifurcation, in a one-dimensional domain with zero flux boundary 
conditions. Fixed parameter values are: Dx = Dy — 2 x 10~ 5 , a = — 1, v = 1, 
and 6=1. The a and /3 parameter values are: a) and d), a = as = — 1, 
/3 = 0; b) and e), a = as = —1.5, /3 = 0.5; c) and f), a — as — —2.0, 
P = 1.0. The total integration time is t = 100, for the time step At = 0.01. The 
one-dimensional region has TV = 2000 lattice sites. In a), b) and c), the initial 
values of X and Y are taken at the unstable focus over the one-dimensional 
lattice except at the center of the domain, where X and Y assume the values 
of the semistable fixed point in phase space. In d), e) and f), the initial values 
of X and Y are taken at the semistable fixed point over the lattice, except at 
the center of the domain, where the system has been perturbed to the unstable 
focus with a finite width perturbation of 160 lattice sites. Propagating fronts, 
traveling waves and aperiodic regimes exist at the SNH bifurcation of the local 
dynamics. 
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Figure 4: Wave length of traveling waves as a function of j3, taken at the SNH 
bifurcation (a — as — f3 — 6), and the same parameter values of Fig. For 
(3 > 0.95, the front develops an aperiodic signal, usually identified as chemical 
turbulence. The variation of the wave length A with f3 follows the approximate 
law, A = c + d/(3. Numerically, d = 0.071 and c = 0.011, for (3 > 0, and 
d = -0.039 and c = 0.024, for (3 < 0. 



for as = —f3 — 1 = —1.5, and we have investigated the limits in a where 
periodic traveling waves persist, Fig. The a parameter controls the distance 
from the stable steady state to the threshold of excitability. After the SNH 
bifurcation, in a range of the a parameter a £ (—1.515,-1.5], traveling wave 
solutions persist, Fig. |5]i). The initial perturbation develops a moving source 
initializing traveling waves, which propagates in the opposite direction of the 
front. When a approaches the value a = —1.515, the amplitude of different 
wave spikes changes, but the wavelength remains constant, Fig. |SJd). Around 
a w —1.515, the propagating signal becomes aperiodic, Fig. [SJ:). Decreasing 
even further a, the turbulent regime behind the moving source decreases. At 
low a values, the moving source is followed by a narrow zone where the generated 
pulse disappears. 

Therefore, after the SNH bifurcation, wave propagating phenomena exists 
in a narrow a parameter range. Geometrically, this parameter range depends 
of the proximity of the saddle and node pair of fixed points. 

For an initial state in the unstable focus and a perturbation at the stable 
fixed point in the phase space, Fig. G3i)-f), depending of the width of the per- 
turbation, the perturbation dies out, or the extended system develops a finite 
train of propagating pulses, or develops an aperiodic regime. In the case of 
the simulations of Fig.[SJi), we obtain a traveling pulse for perturbation radius 
larger than 4 lattice sites. For smaller perturbation sizes, the perturbation dies 
out. For large perturbation widths and decreasing a, the isolated pulses develop 
and propagate in the same region a £ (—1.515, —1.5], Fig. EJi)-e). 

We have analyzed the dependence on the initial perturbation size. Increasing 
the radius of the perturbation, more and more traveling pulses appear, Fig.|S^). 
The perturbation increments, excluding the first one, are equal. For the pa- 
rameter values of Fig. the media can propagate waves with a characteristic 
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Figure 5: Solutions and F(x,t) of the reaction-diffusion equation l|12JI 

after SNH bifurcation, in the same conditions of Fig. In a), b) and c), the 
initial values of X and Y are taken at the unstable focus over all the one- 
dimensional lattice except at the center of the domain, where X and Y assume 
the values of the stable fixed point in phase space. In d), e) and f), the initial 
values of X and Y are taken at the stable fixed point over the lattice except at 
the center of the domain, where the system has been perturbed to the unstable 
focus with a finite width perturbation of 160 lattice sites. Fixed parameter 
values are: D x = D Y = 2 x 1CT 5 , a = -1, v = 1, b = 1.0, (3 = 0.5. The a 
parameter values are: a) and d), a = —1.505; b) and e), a — —1.51; c) and f) 
a = —1.515. For a w —1.515, the generated wave front becomes aperiodic. For 
a < as, waves only persist in a narrow a parameter range. 
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wavelength of 141 lattice sites (Fig.^J, implying that for an increment in the per- 
turbation width of 110% of the wave length, a new propagating pulse appears, 
and we have a moving pulse packet. In Fig. Et>) , we show the time evolution of 
an initial condition leading to two packets of three propagating pulses. For this 
particular choice of initial conditions, if the perturbation is not exactly in the 
unstable focus, then the system can only develop one pulse or the perturbation 
dies out. This situation occurs for perturbation values all over the phase space, 
even using large perturbation radius. 



a) 




T«\ -nil". firn 

Perturbation radius 



b) 





Figure 6: a) Number of traveling pulses as a function of perturbation radius at 
the SNH bifurcation. The initial values of X and Y are taken at the semistable 
fixed point over all the one-dimensional lattice except at the center of the do- 
main, where the system has been perturbed to the unstable focus with a variable 
width perturbation, b) Time evolution of an initial conditions leading to two 
packets of three propagating pulses. In this case the perturbation radius has 250 
lattice sites. All the simulations were carried out in an one-dimension domain 
with N = 4000 lattice sites. The parameters are the same as of Figs. Eb) and 
Et), with (3 — 0.5 and a = as = —1.5. 



The phenomena of persistence of traveling wave trains and pulses at and 
after the SNH bifurcation has a simple explanation. At the SNH bifurcation 
(a = as), when the extended system has initial conditions in the unstable 
focus, any perturbation from it approaches the semistable fixed point. Driven 
by diffusion, the points spatially close to the perturbation generate a line in 
phase space. For (3 ^ 0, this line always intersects the threshold of excitability 
in an infinite number of points and a traveling wave is generated. When the 
extended system has initial conditions in the semistable fixed point, any small 
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fluctuation on the concentrations of X and Y crossing the threshold drives 
nearby points of the extended system to trajectories in phase space with different 
lengths. In this case, we have pulse generation. After the SNH bifurcation, if 
the distance between the stable steady state and the saddle point is small, 
larger perturbations cross the excitability threshold, and an instability in space 
develops. As a decreases away from as, the distance in phase space between 
the stable steady state and the saddle point increases, and the diffusion driven 
instability disappears. 

In two dimensional domains the situation is similar to the one dimensional 
case. For example, preparing the system as in the cases of Figs. |SJa) and (5^), 
we obtain a traveling wave train and a traveling isolated pulse, respectively, 
independently of the parameters of the local dynamics, Fig. |JJi)-b). Further 
numerical simulations for the two-dimensional case shows that the excitable 
local system coupled with diffusion in the two phase space variables can develop 
spirals. However, in the two-dimensional case, after the SNH bifurcation, the a 
parameter range where propagating wave trains and pulses persist differs from 
the one-dimensional case. 




Figure 7: Solutions X(x,t) of the reaction-diffusion equation (|12[1 in a two- 
dimensional domain with zero flux boundary conditions. The local dynamics 
has a saddle node pair of fixed points and an unstable fixed point at the ori- 
gin of coordinates. Blue represents higher concentrations of X and red lower 
concentrations. We also show the spatial variation of the concentration of X 
along a cross section passing through the center of the two-dimensional domain. 
In a), the initial values of X and Y are taken at the unstable focus over all 
the one-dimensional lattice except at the center of the domain, where X and 
Y assume the values of the stable fixed point in phase space. In b), the initial 
values of X and Y are taken at the stable fixed point over the lattice except at 
the center of the domain, where the system has been perturbed to the unstable 
focus with a finite width perturbation of 160 lattice sites. In a), we show the 
formation of a traveling wave. In b), we show a propagating solitary pulse. The 
two-dimensional circular domain is embedded in a two-dimensional square lat- 
tice with 2000 by 2000 sites. The parameter values are: Dx = Dy = 2 x 10~ 5 , 
a = —1, v = 1, b = 1.0, (3 = 0.5 and a = —1.51 < as, as in the cases of Figs.GJ>) 
and 03). In both figures, the total integration time is t = 100. 
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4 Conclusions 



We have obtained a local model of an excitable system exhibiting a Hopf and 
a saddle-node homoclinic (SNH) bifurcation. This model has been constructed 
introducing a phase dependence on the angular coordinate of the versal un- 
folding of the Hopf bifurcation. The bifurcation parameters of the model are 
completely independent. In extended media, the excitability state occurs when 
the SNH bifurcation of the local model is crossed. The excitable local system 
has a hctcroclinic orbit connecting two unstable fixed points, and this orbit is 
close to a stable steady state. 

We have shown that in excitable extended systems with a sharp threshold 
and coupled by diffusion, traveling waves, fronts, isolated pulses, packets of 
pulses and aperiodic wave trains develop. In all the cases, we have taken equal 
diffusion coefficients for the local variables of the model. Depending on the initial 
conditions of the extended system, traveling periodic waves, isolated pulses and 
packets of pulses exist for the same parameter values of the excitable system. 
Propagating pulses and packets of pulses only exist at and beyond the SNH 
bifurcation of the local model. Traveling wave trains exist on both sides of the 
SNH bifurcation, but the wave lengths change discontinuously when the SNH 
bifurcation is crossed. 

Excitable systems with a precise threshold and equal diffusion coefficients 
can show a large variety of behaviors, as found in real chemical systems: Iso- 
lated propagating pulses and fronts, as in the iodate-arsenous acid extended 
system [2], packets of propagating pulses or bursting, as in neurons and pan- 
creatic /3-cells |23, and chemical turbulence [22] ■ Excitable systems with pre- 
cise threshold can also show target patterns similar to those observed in the 
Belousov-Zhabotinskii experiment 2Sj. In the case of the Belousov-Zhabotinskii 
reaction, if we assume that the local kinetics is self-oscillatory and can be de- 
scribed by a two variable model |21| , only large perturbations for the vicinity of 
an unstable steady state can be in the origin of the emergence of target patterns. 
These perturbations can be seen as the local fluctuations of concentration occur- 
ring in the media at the microscopic level. This could explain way pacemakers 
appear sparsely and spontaneously in experiments in the extended media. 

The relation between chemical turbulence associated to the appearance of an 
aperiodic regime as reported here, are dependent of the parameters of the local 
kinetics and are less dependent of the initial distribution of the initial states in 
phase space. 

The model presented here shows qualitatively the known types of topologi- 
cally distinct activity waves as observed in real chemical experiments, and the 
different activity waves can be tuned by different initial conditions, maintaining 
the parameters of the local kinetics fixed. 

Finally, one of the consequences of the approach presented here is that in 
excitable extended systems with sharp threshold (Type II excitability), we can 
have pulse packets with any number of bursts, whereas in models with finite 
width threshold (Type I excitability), packets of propagating pulses have never 
been observed. This enables to design experiments in order to decide between 
the two types of excitability. 
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